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a 

Topological band insulators (TBIs) are bulk insulating materials which feature topologi- 

s 

^Lj cally protected metallic states on their boundary. The existing classification departs from 

O time-reversal symmetry, but the role of the crystal lattice symmetries in the physics of 

these topological states remained elusive. Here we provide the classification of TBIs pro- 

(N 

tected not only by time-reversal, but also by crystalline symmetries. We find three broad 

^ classes of topological states: (a) T states robust against general time-reversal invariant 

(N 

perturbations; (b) Translationally-active states protected from elastic scattering, but sus- 
^ ceptible to topological crystalline disorder; (c) Valley topological insulators sensitive to 

the effects of non-topological and crystalline disorder. These three classes give rise to 18 

•i-H 

different two-dimensional, and, at least 70 three-dimensional TBIs, opening up a route for 



the systematic search for new types of TBIs. 
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Topological phases of free fermionic matter are in general characterized by an insulating 
gap in the bulk and protected gapless modes on the boundary of the system[[Ill3. Integer quan- 
tum Hall states represent first examples of topologically protected phases in absence of any 
symmetries with the topological invariant directly related to the measured Hall conductance^. 
Recently, it became understood that even in the presence of fundamental symmetries such as 
time-reversal, topologically protected states of matter can, in principle, exist. In particular, it 
has been shown that time-reversal invariant (TRI) insulators in two dimensions (2D)(4l and 
three dimensions (3D) [|51 |6l are characterized by Z 2 topological invariants which pertain 
to the existence of the gapless boundary modes robust against time-reversal preserving per- 
turbations, and may host Majorana quasiparticles[8J, as well as provide the condensed-matter 
realization of the theta-vacuum[|9]|. The theoretical prediction [fTUl ITTil and experimental real- 
ization of the Z 2 -invariant topological band insulators lfT2l IT3l [T4l [T5l [T6l gave a crucial boost 
in the understanding of these exotic phases of matter which culminated in the general classi- 
fication of topological insulators and superconductors based on time-reversal symmetry (TRS) 
and particle-hole symmetry (PHS) within the so-called tenfold periodic table llT7l[T8l[T9l . 

The role of the crystal lattice in this classification is to provide a unit cell in the momen- 
tum space, the Brillouin zone (BZ), topologically equivalent to the <i-dimensional torus, over 
which the electronic Bloch wavefunctions are defined. The tenfold classification follows then 
assuming that all the unitary symmetries of the corresponding Bloch Hamiltonian have been 
exhausted and therefore the only remaining symmetries are, according to the Wigner's theo- 
rem, antiunitary TRS and PHS. In three dimensions (3D), however, by considering a Z 2 TBI 
as a stack of two-dimensional (2D) ones, thus assuming a layered 3D lattice, additional three 
"weak" invariants associated with the discrete translation symmetry have been foundjH, which 
characterize these states, and may be probed by lattice dislocations [|20l . On the other hand, it 
has been realized that topological states protected just by lattice symmetries, such as inversion, 
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can exist in 3D EUE21 1231 EH 1221 ISO- Additionally, it has been recently found that a state pro- 
tected both by TRS and the lattice C4 rotational symmetry susceptible to the lattice dislocations 
can be possible in 2D[27J. However, in general, the role of the space group symmetries in the 
physics of topological states remained elusive. 

We here provide the complete classification of TBIs protected not only by TRS, but also by 
space group lattice symmetries. As a starting point, we depart from the construction by Fu and 
Kane [0171 ETJ to compute the Z 2 invariant in terms of the matrix of overlaps 

Wmn = (Wm(-k)|$K(k)}, (1) 

where ■d is the time -reversal operator and |w n (k)) is the n-th occupied Bloch wavefunction. The 
quantities of central significance are 

_ VdetHrp] 

1 PfKr,)] {) 

defined at the points Tj in the BZ where the Hamiltonian commutes with the time-reversal 
operator. Since the matrix w is antisymmetric at the points Tj, the Pfaffian is defined at these 
points and detfuifT;)] = (Pf[w(r,)]) 2 . The topological Z 2 invariant, v, is then given by (—1)" = 
Yi T , Si, and its non-triviality implies a topological obstruction for defining the wavefunctions 
through the entire BZ with an even number of band inversions. Notice that the evaluation of the 
topological invariant in terms of the signs of the Pfaffian does not depend on the dimensionality 
but only on the fact that the Hamiltonian possesses TRS which, due to the vanishing of the 
Chern number, guarantees the existence of globally defined wavefunctions throughout the BZ. 

First, notice that the set of the points Ti at which the Hamiltonian commutes with the time- 
reversal operator is fixed by the space group of the lattice, see Table [TJ Second, we choose the 
overall phase of the Bloch wavefunctions so that a unique phase, which we dub the 'T" phase, 
has Sr = —1 at the T-point in the BZ and Si — 1 at all the other high symmetry points. A 
crucial observation is that the distribution of signs of the Pfaffian, <5,, at the points Tj, and not 
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only their product, encodes the additional topological structure. To show this, we first consider 
how the matrix of overlaps transforms under a lattice symmetry operation represented by a 
unitary operator U 

w mn (k) = (u m (-k)|% n (k)> = (u m (-Uk)\mU^\u n (Uk)) = w mn (Uk). (3) 

As a consequence, when some of these high symmetry points are related by point-group sym- 
metry of the lattice, their signs of the Pfaffian have to be equal. Therefore, it is sufficient to 
consider a subset, T a , of representative, inequivalent high symmetry points that are also not 
related by any symmetry. This leads to the following rule that allows for determination of all 
the topological phases given the space group and the corresponding high symmetry points, Tf. 
each phase is obtained by selecting a single representative high-symmetry point T a and set- 
ting Sr a = —1, which automatically sets 5r b = —1 at all the high-symmetry points Tb related 
by point group symmetry to T a . Such phases are separated by a topological quantum phase 
transition that involves bulk bandgap closing which changes the values of 5j's. 

Let us illustrate this simple classification principle by an elementary example. We start with 
the T phase on a square lattice, 5 r = —1, and 5 X = 8y = &M — 1 where X, Y, and M are 
the time-reversal invariant (TRI) momenta in the BZ. By applying our rule, we immediately see 
that, besides the T phase, we obtain an "M" phase with 5m = — 1, and 5r = Sx = 5y = 1 (Table 
[T]). This phase is disconnected from the T phase through a topological quantum phase transition 
with the bandgap closing at the X and the Y points. This phase is protected by TRS and 
is also susceptible to dislocations ffTTTl . and represents an example of a "translationally-active" 
phase. Furthermore, since the X and the Y points are related by a C 4 rotation, there can exist 
a phase with Sx = Sy = — 1, and Sr = 5m = 1- The product of the S^s at all TRI momenta 
then yields the trivial Z 2 invariant, v = 0. However, C 4 rotational symmetry protects this 
phase, since it pins the band inversions at the X and Y points. This phase represents a "valley" 
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(or "crystalline" E3l ) insulator - a phase trivial tenfold way-wise but protected by the lattice 
symmetries. Indeed, this new phase, which we dub the "X-Y" phase, is realized in an extension 
of the M-B model for a quantum spin Hall insulator with the next-nearest neighbor hopping, as 
shown in Fig. [T| and responds nontrivially to dislocations; see Supplementary Information for 
the technical details. When C 4 rotational symmetry is reduced to C 2 , the X and Y points are no 
longer related by symmetry, and therefore the symmetry constraint on the 5x and Sy is no longer 
present. We then expect the X-Y phase to be unstable, and to yield instead nontrivial phases 
with 5 X = — 1 or 5 Y = — 1, and Si — 1 at all other TRI momenta. Our calculations, indeed, 
confirm this within the M-B tight-binding model, as shown in Supplemetary Information A. In 
general, an even number of TRI momenta related by symmetry yield a valley phase, protected 
by crystal symmetry while having v = 0. 

Let us now elaborate on the role of the space group of the underlying lattice in this clas- 
sification, since this symmetry group defines the relation between the high symmetry points. 
The difference in phases found on rectangular and rhombic lattices serves as a clear illustration. 
Both these lattices have D 2 point-group symmetry, but different wallpaper groups (space groups 
in 2D). The rhombic case has two inequivalent TRI momenta related by point group symmetry 
and hence a valley phase, see Table [T] and Supplementary Information, Sec. A. On the other 
hand, in the rectangular case all D 2 symmetry operations map any TRI momentum to its equiv- 
alent, thus no valley phase is possible. From Table [T] it is seen that in 2D the phases, as related 
to space groups, cluster in Bravais lattice classes, with an exception: the hexagonal lattice. We 
will see that this clustering is less generic in 3D. In turn, the primitive Bravais hexagonal (tri- 
angular) lattice ipQmm) is invariant under the C 6 rotational symmetry around a lattice site, as 
opposed to the non-primitive hexagonal lattice (p3ml) realized in graphene. Based on our rule, 
we conclude that in the latter case only the Y phase is possible, which is in fact realized in the 
Kane-Mele model ll28ll . In contrast, on the former lattice (triangular), the points K + and K_ 
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are related by a C§ symmetry and thus each of these points becomes TRI. The number of TRI 
momenta is increased, ultimately yielding a possibility of additional translationally-active and 
valley phases, as shown in Table [T] These phases are realized within the M-B tight-binding 
model (Supplementary Information, Sec. A), and their robustness against disorder is shown in 
Supplementary Information, Sec. B. 

The above rule allows us to completely classify and index the topological phases: the last 
entry in Table [T] The set of BZ high-symmetry points I" 1 , at which there is band inversion, 
i.e., 8 Vi = —1, is invariant under the operations of a subgroup of the lattice space group. This 
symmetry subgroup therefore protects and labels the topological phase. The other element in 
this indexing is the protection by TRS (T), existing when Z 2 invariant v — 1, giving, for instance 
T — pAmm as the Y phase on the square lattice. When the protecting symmetries coincide 
between phases, we explicitly label Tj (lower index), as, e.g., for T — p2mx, T — p2my and 
T — p2m M phases on the rectangular lattice. This leads to the list of topological phases in 2D 
presented in Table[T] which gives 18 distinct topological phases. As our general result, there are 
two additional broad classes of topological states protected by TRS or crystalline symmetries, 
besides the class of states robust against general TRS perturbations (r-states): translationally- 
active states protected both by TRS and lattice symmetry, responding to dislocations, and valley 
insulators which are tenfold- way- wise trivial but protected by space group symmetry and also 
susceptible to dislocations. 

Our procedure can be applied in the same way in 3D but it becomes more involved given the 
230 space groups and the large number of high-symmetry points. We find at least 70 different 
phases OUl . Here we will illustrate these matters for a number of simple crystal structures 
(Table[2]) which include those of TBIs that are of present empirical relevance |[P3l IT4l [T5l [T6ll . To 
illustrate matters, consider the primitive cubic lattice (Table[2]) with the familiar eight TRI points 
(Fig. [2]\). Crucially, the points (X,Y, Z) are related by a three-fold rotation, as well as the 
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points (X', Y', M). Consequently, we obtain four TRS protected phases. We notice that this is 
quite different from the indexing procedure introduced by Moore and Balents [5J. For instance, 
our T — pm3m (T) and T — p3(4) R (R) phases correspond with their (1; 0, 0, 0) and (1; 1, 1, 1) 
indices, respectively. Their latter two indices would also correspond with the T — £>3(4)m 
and T — p3(4) x phases, respectively. The other possibilities in their classification are either 
coincident with our four TBIs, or represent a 3D phase not protected by crystal symmetries due 
to implicit dimensional reduction (e.g. layered 3D lattice); see Figs.[2]A.,B. 

The power of the space group classification becomes further manifest for non-cubic lattices. 
Consider the 3D hexagonal lattice which consists of two hexagonal layers with the wallpaper 
group p6mm stacked on top of each other. The TRI momenta comprise two copies of the ones 
on the 2D hexagonal lattice, separated by a perpendicular translation. Accordingly, the phases 
can easily be obtained by considering the k z = plane (Table [T]), which contains the V, M 
and K points, and those of the other translated plane associated with the points A, L and H, 
respectively (Table [2]). Consequently, there are eight TRS protected phases resulting from the 
combinations of a TRS protected phase in one plane and a trivial or a valley configuration in the 
other plane. Additionally, there are two valley phases which are configurations with one plane 
featuring a valley phase and the other a trivial configuration. Notice that a potential double 
valley phase with a valley phase in each of the planes is not protected by 3D crystal symmetry 
and is therefore trivial. We again point out that the truly 3D valley phase is determined by a 
three-dimensional point group, i.e., the one whose action cannot be reduced to the 2D case. 

Most experimentally observed TBIs are of the T kind, such as Bi 2 Se 3 [|T5l and Bi 2 Te 3 [[T6|, 
except for Bi z Sbi^ where the 3D pursuit started |[T3l [141 which is r3m — T — r3 L . The 
theoretically predicted rocksalt actinides 11291 are actually of the translationally-active class 
fm3m — T — / 3(4). SnTe, as well as Sn-doped compounds PbTe and PbSe have the same space 
group but the phase recently observed in these compounds has no TRS protection OTl l32l l33l 
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and is a valley phase [|23l . Let us inspect this phase in more detail. This phase turns out to be 
indexed as fm3m — /3(4) (Table [I]). The system has mirror planes in the momentum space 
formed by the T and any two of the L points, which thereby relate the remaining two L points 
by symmetry. As a result, a mirror-symmetric crystal cut along T — X — T line in the (001) 
surface features a pair of Dirac cones (a double Dirac cone) which is therefore also protected 
by the same symmetry ll24ll . Notice that in addition we predict valley phases at the W and the 
U points in the Brillouin zone protected by both the four-fold and the three-fold rotational sym- 
metries, labelled by fm3m — fA3 w and fm3m — fA3u, respectively. The P^-phase originates 
from six inequivalent symmetry-related W -points in the BZ where a band inversion gives rise 
to a valley phase. In addition, in the same phase, the (001) surface features Dirac cones for the 
cut along r — M — f and f — S — f lines, but not along f — X — f direction, as it is the case 
in the fm3m — /3(4) phase, see Fig. 2C. Therefore, the detection of the Dirac cones in the 
f — M — f and f — S — f directions in ARPES experiments would be a clear signature of this 
valley phase. 

In conclusion, we provided the space group classification of the TBIs both in 2D and 3D. 
As a result, we found two additional broad classes of topological phases, besides the TRS pro- 
tected T states: the translationally-active phases, protected by both TRS and crystal symmetries, 
but susceptible to topological crystalline disorder, and the valley phases solely protected by the 
space group symmetry, and therefore susceptible to both elastic and topological crystalline dis- 
order. Our complete classification scheme based on the full 2D and 3D space groups has as most 
important consequence that it demonstrates the potential existence of at least seventy distinct 
topological phases of insulating matter and we anticipate that this will be a valuable guide in 
the future exploration of this landscape. 
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1 Tables 



Bravais Lattice (PG) 


WpG 




Si 


Index (Phase) 


Square (D 4 ) 


pAmm 
pAgm 
pA 


(r, M, X, Y) 


(-1,1,1,1) 


T-pAmm (F) 


(1,-1,1,1) 


T-pA (M) 


(1,1,-1,-1) 


pA (Jf-r-valley) 


Rectangular (D 2 ) 


p2mm 
p2mg 

pZgg 

pm, pg 


(Y, M, X, Y) 


(-1,1,1,1) 


T-p2mm (F) 


(1,-1,1,1) 


T-p2m M (M) 


(1,1,-1,1) 


T-p2m x (X) 


(1,1,1,-1) 


T-p2m, Y (Y) 


Rhombic (D 2 ) 


c2mm 
cm 


(r, M , M_i, Mi) 


(-1,1,1,1) 


T-c2mm (T) 


(1,-1,1,1) 


T-c2m, (M ) 


(1,1,-1,-1) 


c2m (M-valley) 


Oblique (C* 2 ) 


P 2 
pi 


(r, M , M_i, Mi) 


(-1,1,1,1) 


T-p2 (D 


(1,-1,1,1) 


T-p2 Mo (M ) 


(1,1,-1,1) 


T-p2 ALl (M_i) 


(1,1,1,-1) 


T-p2 Ml (Mi) 


Hexagonal 

(hexagonal - Dq) 

Hexagonal 
( rhombohedral - D3) 


pQmm 
p6 


(r, M , M_ l5 Mi, K_, K+) 


(-1,1,1,1,1,1) 


T-pQmm (T) 


(1,-1,-1,-1,1,1) 


T-p6 (M) 


(1,1,1,1,-1,-1) 


p6 (if-valley) 


p3ml 
p3lm, p?> 


(r, M , M_i, M x ) 


(-1,1,1,1) 


T-p3ml (T) 



Table 1: Table of the topological phases in two dimensions. For each of the lattice structures, the 
corresponding point-group (PG) symmetry and the relevant wallpaper group (WpG), i.e. space 
group, are given. The corners of the square and rectangle are denoted by M, whereas in the 
triangular Bravais structure they are indicated by K. Additionally, the centers of the edges are 
denoted by X and Y in both the square and rectangular case and by M in the other lattices [34]. 
The resulting phases are characterized by the distribution of 5i at the Tj points consistent with 
the WpG symmetry. Phases cluster in Bravais lattices, with the hexagonal structure being the 
only exception. In this case the WpGs containing six-fold and three-fold rotational symmetries 
relate the high symmetry points in different ways. As a result, the Hamiltonian does not com- 
mute with the time-reversal operator at the K points in the latter case. The obtained phases are 
ultimately protected by TRS (whenever v = 1), WpG symmetry, or both, and are accordingly 
indexed. The index (last column) describes the part of wallpaper group that leaves the subset T, 
having 5i — — 1 invariant, while the additional label 'T' denotes TRS protection. In the column 
denoted "Phase" we introduce a convenient but imprecise shorthand notation. 
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Bravais Lattice 


PGS 


SG 




Si 


Index (Phase) 


Primitive Cubic 


o h 


pm3m 
pn3n 
pn3m 
pm3n 


(r, R, X, M) 


(-1,1,1,1) 


T-pm3m (T) 


(1,-1,14) 


T-p3(4) R (R) 


(1,1,-1,1) 


T-p3(4) x (XYZ) 


(1,1,1,-D 


T-p3(4) M {MX'Y') 


Hexagonal 




p6mm 
p6cc 
pQ^cm 
pQ^mc 


(T, M, A, L, K, H) 


(-1,1,1,1,1,1) 


T-p6mm (T) 


(1,-1,1,1,1,1) 


T-p6 M (M) 


(1,1,-1,1,1,1) 


T-p&a (A) 


(1,1,1,-1,1,1) 


T-p&l (L) 


(-1,1,1,1,1,-1) 


T-p&th (TH) 


(1,-1,1,1,1-1) 


T-p& H M (MH) 


(1,1,-1,1,-1,1) 


T-p& KA (KA) 


(1,1,1,-1,-1,1) 


T-p6 LK (LK) 


(1,1,1,1,-1,1) 


p6x (if-valley) 


(1,1,1,1,1-1) 


p6H (i?-valley) 


Face Centered Cubic 


o h 


fm3m 
fm3c 
fd3m 
fdlc 


{V, X, L, U, W) 


(-1,1,1,1,1) 


T-fm3m (T) 


(1,-1,1,1,1) 


T-/3(4) (X) 


(1,1,-1,1,1) 


/3(4) (L-valley) 


(1,1,1,-1,1) 


/43c/ ([/-valley) 


(1,1,1,1,-D 


f43 w (W-valley) 


Rhombohedral 


D 3d 


r3m 
r3c 


<T, L, F, Z, P, K, B) 


(-1,1,1,1,1,1,1) 


T-r3m (T) 


(1,-1,1,1,1,1,1) 


T-r 3 L (L) 


(1,1,-1,1,1,1,1) 


T-r3 F (F) 


(1,1,1,-1,1,1,1) 


T-r3 z (Z) 


(1,1,1,1,-1,1,1) 


r3p (P-valley) 


(1,1,1,1,1,-1,1) 


t3k (if -valley) 


(1,1,1,1,1,1,-1) 


r3 B (S-valley) 



Table 2: Topological phases anticipated in 3D for some specific point group (PG) symmetries. 
Bravais lattices with same PG symmetries have different space groups (SG). We point out that, 
in contrast to the 2D case, the phases do not cluster in Bravais lattice structures. For example, 
the fourfold rotational symmetry crucial for the /43[/ and /43yj/ phases is not contained in every 
space group associated with the face centered cubic lattice. As cVs attain the same value at the 
points related by lattice symmetry or a reciprocal lattice vector, only one representative is 
given from each set of such points. We note that the rhombohedral T-t3l phase is observed in 
BuSbi_^ lfi~3l [i"4l . while the T-r3m phase is found in Bi 2 Se 3 lfi~5l and Bi 2 Te 3 [[16]|. Moreover, 
the fm3m — /3(4) phase has recently been observed in SnTe [f311 , as well as in Sn-doped 
compounds PbTe [32| and PbSe (33). 
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2 Figures 




Figure 1 : Phase diagram of the extended M — B tight-binding model. As function of the model 
parameters M/B and z — B/B, where B (B) is the (next) nearest-neighbor hopping parameter 
and M the difference in on-site energies, the different distributions of Si are obtained with the 
corresponding phases listed in Table [TJ consult SOM A for details. Furthermore, the spectra 
of edge states per spin component are displayed for the non-trivial phases, demonstrating that 
the valley X-Y phase exhibits a pair of Kramers pairs of metallic edge states. The real space 
localization of these edge states is also presented, where the radii of the circles represent the 
magnitude of the wave function and the colors indicate the phases as shown on the left. 
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Figure 2: Illustration of the role of lattice symmetries in the classification of topological states. 
(A) The eight TRI momenta in the Brillouin zone of the primitive cubic lattice. When only TRS 
is considered the sign of any quadruple of <Vs within a plane connecting them can be changed, 
leaving their product the same. As a result one obtains, in addition to the 'strong' invariant, three 
weak invariants corresponding to the orthogonal planes. (B) The constraints on the <Vs arising 
from the lattice symmetries. The high symmetry axes A, A and £ represent axes of four-, three- 
and two-fold rotations, respectively; these transform the TRI points in the colored planes into 
each other, and thus constrain the corresponding 5;'s to be equal. (C) The Brillouin zone of the 
face centered cubic lattice with high symmetry points and a mirror plane that projects onto the 
r — X — T line in the (001) plane. The VT-valley phase features Dirac cones along f — M — f 
and r — S — T lines, but not along F — X — t lines. 
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A. Model details 

We present details on the tight-binding models used to illustrate the phases in two dimensions. 
Specifically, we consider models of the form as introduced in Ref. [SI] which we generalize to 
the five Bravais structures in two dimensions. These models assume nearest neighbor hopping 
on the specific lattice, which consists of spin degenerate \s) and \p x + ip y ) type orbitals on every 
single site, and are therefore, in natural units U — c — e — 1, described by 



Here {6} denotes the vectors connecting the nearest neighbors on a particular Bravais lattice 
and ^ = (s t (R), p t (R), s^R^p^R)) annihilates the s and p type orbitals at site R. We set 
the hopping matrices for the spin up and spin down electrons 



with tg = — 1A5 exp(i(pg) and t' s = —iAsexp(—i(ps), in terms of the polar angle ipg of the 
vector 8, A s / P)(5 = D ± B§, and on-site energies e = [(C — 4_D)r + (M — 4£?)r 3 ] <g> a , where 
the Pauli matrices r and cr act in the orbital and spin space, respectively. Here, t and cr are 
the 2 x 2 unity matrices. 




(SI) 
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By performing a Fourier transform, the Hamiltonian dST) assumes the form 



«TB = £*t("i k) „»,„],„ ,S2, 

k 



H*(-k) 

K 

where 

i7(k) = r • d(k) 

and the vector d(k) obviously depends on the form of the lattice. The upper and the lower block 
of this Hamiltonian are related by the time-reversal symmetry represented by an antiunitary 
operator T = r <g> ia 2 K, with K as the complex conjugation. The form of the tight-binding 
Hamiltonian in momentum space is quite generic in a sense that it is valid also when longer- 



range hopping terms, such as next-nearest neighbor ones, are included in the Hamiltonian (SI ), 
as we will see later. As a result, the spectrum of the above Hamiltonian takes the simple form in 
terms of the vector d, E(k) = ±V d ■ d = ±|d|. The closing of the bandgap occurs when the 
magnitude of the vector d(k) vanishes, i.e., when d(k) = 0, which is the case only at the special 
high-symmetry points in the Brillouin zone (BZ) and for special values of the parameters in the 
Hamiltonian, as we will see in the specific examples. 

Taking up the square lattice and neglecting the energy shifts, we obtain the familiar three 
component d(k) vector with 0^,2 (k) = ± sin(/c Xi j / ) and d 3 = M — 2B{2 — cos(k x ) — cos(k y )). 
The corresponding spectrum, doubly degenerate in the spin space, is gapped except for the 
values M /B = 0, 4 or 8. When < M/B < 4 the band gap opens at the T point in the BZ and 
the system enters a topological phase, denoted as the T phase. For 4 < Mj B < 8 the system is 
still topologically non-trivial (M phase), but with the band gap opening at k = (it, 0) (X point) 
and k = (0, tt) (Y point). In contrast, the other values of M/B give rise to a topologically 
trivial phase. Interestingly, the resulting nontrivial phases in two dimensions can (if needed by 
adiabatic continuity) be visualized by the position of the Skyrmion of the band structure field 
d(k) = d(k)/|d(k)| with the associated density s(k) = d(k) • [^^(k) x dk y d(k)}. Moreover, 
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these phases can be probed by dislocations, and only the M -phase results in zero-modes bound 
to the dislocations, whereas a 7r flux binds zero-energy modes in both non-trivial states [52]. 

The connection with the characterization presented in the main text is achieved by relating 
the Pfaffian to the band- structure vector d(k) through 6i = sign[d 3 (Ti)], which can be readily 
shown, as well as by noting that Tj comprise the points at which the gap closes. As a result 
the T and M phase are indeed characterized by an opposite sign of the Pfaffian at those points. 
Furthermore, the C4 rotation connects the 5j's and thus dictates that the values at the X and Y 
points have to be the same. As a result, the 5i's can assume a configuration for which 5 X = 
Sy 7^ 5r = 5m, which we call the X — Y valley phase. This phase can easily be incorporated in 



the model ( [ST] ) by introducing a next nearest neighbor term with the parameters defined exactly 



as in the nearest neighbor case. Accordingly, we obtain 

(sm(k x ) + cos(k x ) sin(k y ) 
— sm(ky) — sin^.) cos(k y ) | . (S3) 

M - 2B[2 - cos(k x ) - cos(k y )} - 4B[1 - cos{k x ) cos{k y )} / 

The gap at the TRI momenta and hence the configuration of 5i can be tuned as function of the 
parameters B/B and M/B resulting in the phase diagram shown in Fig. 1 in the main text. 
Importantly, we find that the resulting valley phase features two pairs of edge states (Fig. 1 
in the main text). In this case there are two single Skyrmions located at the X and Y points 
and integration of the Skyrmion density yields a winding number two per spin. Moreover, this 
phase can also be probed in the exact same manner as the T and M phases (see Supplementary 
Information, Section B). 

As the next step, the rectangular model is obtained by making the magnitude of the hop- 
ping parameters anisotropic, which effectively reduces the fourfold rotational symmetry to the 
twofold, see also Section B in Supplementary Information, Eq. ( |S8| ). Consequently, the oppo- 
sitely valued Si's can be located at any of the TRI momenta T, M,X, Y, and the corresponding 
phases can be distinguished by insertion of dislocations (see Supplementary Information, Sec- 
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Figure SI: The Brillouin zone of the primitive hexagonal (triangular) lattice and the phase 
diagram of the M-B tight-binding model on the same lattice. High-symmetry points in the 
Brillouin zone at which the Hamiltonian commutes with the time-reversal operator are indi- 
cated. Additionally, the Chern number per spin block (C = = —C^) is shown for each 
phase. 

tion B). We find that by increasing the ratio r of hopping magnitudes in the y and x directions 
the valley phase is lost: the system either enters a trivial state, a translationally-active X phase 
or a phase with C = 2 per spin not protected by either TRS or lattice symmetry. Namely, only 
the former two phases are realized when the ratio r is large enough. The two edge modes per 
spin in the X-Y valley phase are gapped out upon an infinitesimal breaking of C4 symmetry, 
i.e. by r > 1. The phase diagram of the model for r = 2 is shown in Fig. S5. 

In case of a primitive hexagonal (triangular) lattice with the wallpaper group p6mm and 
with the Brillouin zone shown in Fig. SI, the nearest neighbor hopping part of the Hamiltonian 
is described by 

\ 



d(k) 



( sinfca. + 2sin(^) cos(7r/3) caa{&k y ) 
-2 cos(^) sm(vr/3) sin{^k y ) 



\M - 2B[2 - cosk x - 2 



cos 



cos 



^3 



(S4) 



k y )\.) 



As for the other Bravais lattices, the sign of the Pfaffian at the Fj points is Si = signfc^^)] and 
we observe that due to the very presence of the lattice symmetry (C 6 in particular), Si is also 
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Figure S2: The skyrmion lattices in the extended Brillouin zone for the three topological phases 
on the primitive hexagonal (triangular) lattice. The locus of the Skyrmions is colored in red. 



(a) The Skyrmion lattice for the T phases within the model ((ST), (b) K valley phase within the 
same model, (c) When next nearest neighbor hopping, z = B/B, is taken into account, the 
system can enter the M phase, which has the Skyrmions positioned at the M points. 

well defined at the two inequivalent corners K + and K- of the Brillouin zone. Additionally, the 
symmetry relates the centers of the edges M and we thus expect the phases as indicated in the 



Table 1 in the main text. The resulting Hamiltonian corresponding to the vector d in (|S4]) indeed 
exhibits the T and the K valley phase for —2 < MjB < 6 and 6 < M / B < 7, respectively. 
The M phase can also be captured by including the next nearest neighbor hopping term 



Asm(^k x ) cos(7r/6) cos(^-k y ) ^ 
ci(k) - 1 -A[sm(y/3k y ) -2 cos^) sin(7r/6) sin(^fc„)] 
-2B[3 - cos(VSk y ) -2cos(^k x ) cos(^k y )] j 



(S5) 



For strong enough next-nearest neighbor hopping, z = B/B > |, the gap closes at the K points 
before closing at the M points, for increasing M/B, and as a result the system exhibits a T and 
M phase, see Figs. SI and S2. We note that we need six inequivalent high symmetry points, 
marked in Fig. SI, to capture the three observed phases within the rule presented in the main 
text. In contrast, the symmetry relating K + and K_ can readily be broken by considering a 
non-primitive honeycomb lattice structure, realized in graphene, which then makes K± related 
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by TRS. The Hamiltonian does not commute with the time-reversal operator at these points, 
and therefore, these two points no longer belong to the set IV As the M points are still related 
by symmetry, we expect only a T phase in this case. A prominent example is the Kane-Mele 
model [S3] with the configuration of Si's explicitly calculated in Ref. [S4] using a rhombic unit 
cell. Drawing this cell in the extended BZ we observe that this configuration indeed corresponds 
with a T phase in the usual Wigner-Seitz cell. Additionally, the result can also be understood 
from the Skyrmion picture. Due to inversion-asymmetry the gap closes at K± points instead of 
closing at the TRI momenta [S5, S6]. The spin sub-blocks exhibit merons at K + and K_ points 
which can be explicitly shown by rewriting the Pontryagin index as in Ref. [57]. 

Finally, we consider the oblique and rhombic Bravais structures distinguished by the pres- 
ence of mirror symmetry in the latter case. This symmetry is ultimately responsible for the 
existence of the valley phase on the rhombic lattice. Considering the rhombic lattice, it can be 
readily shown that 

(A\ cos(0i) sin(^) cos(k y ) + A2 sin(rik x ) \ 
—A\ sin(0i) sm(A; x ) cos(k y ) — A3 sin^fcy) , 
M - 4B[1 - cos{k x ) cos(^)] + 2B[t{ 1 + tvT 1 - rf 1 cos^i/^) - tvT 1 cos{r 2 k y )}) 

(S6) 

where k x = e\ cos((pi)k x , k y = e\ sin((/>i)/c y , with 0i as the polar angle of the vector ei, and 
ei = |ei I ; see Fig. S3. Furthermore, A { is the hopping between the s and the p orbitals along the 
vector while B is the hopping between the two s orbitals along the vector ei, and the hopping 
between these orbitals along the vectors e 2 ,3 is -B 2 ,3 = Bt^. Due to the mirror symmetry, the 
M_i and M 1 points are connected and as a result the gap closes at these points for the same 
value of M/B. Consequently, tuning the gaps at the TRI momenta by varying the parameters 
of the model, the configurations of ^'s and the corresponding phases shown in Table 1 in the 
main text are easily obtained. The oblique lattice can be then treated similarly. Denoting the 
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Figure S3: The left panel shows the real space rhombic lattice with the vectors connecting the 
neighbors. The right panel shows the corresponding Brillouin zone with the TRI momenta for 
e2 = | e 2 1 =3 and 63 = |e 3 | = 4 in arbitrary units. Crucially, the mirror symmetry connects the 
M_i and Mi points. 

different connections between neighbors as indicated in Fig. S4, we obtain 

(Ai cos((f>i) sin(ei • k) + A 2 cos(</> 2 ) sin(e 2 • k) + A 3 cos(4>3) sin(e 3 • k) + A 4 cos(</>4) sin(e4 • k)\ 
— A\ sin(0i) sin(ei • k) — A 2 sin(^ 2 ) sin(e 2 • k) — A 4 sin(</>4) sin(e4 • k) 
M — 2B[2 — t^ 1 cos(e! ■ k) — r 2 _1 cos(e 2 • k) — t^ 1 cos(e 3 • k) — t^ 1 cos(e 4 • k)] J 

(S7) 

where angles i5 i = 1, ...,4, are the polar angles of vectors ej. The s — p hopping along 
the vector is Ai, and the s — s hopping along the same vector is = Bt[~ ; see also Eq. 
(SI). Again, by varying the model parameters (consequently also the shape of the BZ), it is 
straightforward to realize the different phases which are also visualized by the presence of a 
Skyrmion at either of the indicated points. More interestingly, the analysis confirms yet again 
the significance of the symmetry group of the underlying crystal. Essentially, the only difference 
between the oblique and rhombic systems is the presence of the mirror symmetry. Therefore, 
the same (1, 1,-1,-1) configuration of Si's results in a valley phase with the Chern number 
equal two per spin block only when the mirror symmetry is present, i.e., only in the case of the 
rhombic lattice. 
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M ■ 



Figure S4: The real space oblique lattice and its corresponding Brillouin zone. In contrast to the 
rhombic case, the TRI momenta are not related by a lattice symmetry. In this example, vectors 
ei and e 2 connect nearest-neighboring sites, while e 3 and e 4 connect next-nearest neighboring 
sites. 

B. Probing the topologically non- trivial phases 



As already shown by the present authors within the M-B tight-binding model on the square 
lattice, a 7r-flux probes Z 2 topological order through binding of the zero modes, while the dis- 
locations act as universal probes of the translationally-active topological states [52]. Here we 
show that indeed the same holds for the valley phase on the square lattice, realized within an 
extended M-B tight-binding model which includes next-nearest neighbor hopping. Moreover, 
we demonstrate by considering the same model on the rectangular lattice that the directionality 
of the dislocation can serve to distinguish translationally-active topological phases at different 
high symmetry points in the BZ. 



The M-B tight-binding model ( [ST] ) on a rectangular lattice assumes the form (S2) in terms 
of the vector 



/ 



d(k) 



svn{k x ) + ^7j=p^ cos(<£v)[sin(A; + ) — sin(fe_)] 
— 1/r sm(rk y ) — 77^= sm(ip r )[sm(k + ) + sin(A;_)] 



\M - 2B[2 - cos(k x ) - 1/r cos(rk y )\ - 2B[2 - cos(k + ) - cos(>_)]/ 



(S8) 



where r denotes the ratio of the lattice spacing in the y and the ^-direction, ip r = arctan(r), 
r = A x /A y = B x /B y , and k± = ±k x + rk y . In the remainder we conveniently fix r = r = 2 
and A — 1. The resulting phase diagram is presented in Fig. S5. 

Considering first only nearest neighbor hopping, the non-trivial phases on both the square 
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and rectangular lattices include the T and M phases. Accordingly, we find that the M phase 
gives rise to a dislocation zero-mode, whereas a n flux results in a zero-mode in both the V and 
the M phase (Fig. S6). The dislocation zero mode hybridizes with the low-energy edge states 
in a finite system. Separating these states by introducing weak disorder via a random chemical 
potential, we confirm that the edge state indeed is at a finite momentum k = ir. 

The next nearest neighbor hopping, as already shown, realizes the valley phase on the square 
lattice (Fig. 1 in the main text). As anticipated, introduction of a magnetic n flux results in a 
pair of Kramers pairs of zero modes in the system. More interestingly, the bulk-boundary 
correspondence and the directional character of a dislocation result in the possibility of dis- 
criminating between the two valley points. Namely, a dislocation with Burgers vector b along 
the x or y direction binds only a single pair of zero-energy modes. Analysis of the edge state 
hybridized with the dislocation zero-mode then confirms that the dislocation indeed "sees" the 
valley at the momentum parallel to the direction of the Burgers vector (Fig. S7). 

When r is increased (breaking square symmetry to rectangle), the Chern number per spin 
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Figure S5: Phase diagram of the extended M — B model on the rectangular lattice in Eq. ( |S8] ) 
for r = r = 2 and A = 1. The nontrivial phases include the T, M and a X phases. Additionally, 
the Chern character C per spin, C = Cf = — Q. is indicated. 
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Figure S6: Effect of introducing a n flux or dislocations in the T and M phase of the square lat- 
tice model. The top row shows two zero-modes bound to 7r-anti-7r flux pair (periodic boundary 
conditions) appearing in the spectrum in case of the T phase and the real space localization of 
these modes. For convenience we only consider one spin component. The bottom row shows the 
effect of a dislocation-anti-dislocation pair (periodic boundary conditions) with Burgers vector 
b = e y in the M phase. The lattices contain 33 x 33 unit cells. 



block may reduce to zero and insertion of a 7r-flux or a dislocation no longer results in a zero- 
mode in the system, confirming the intimate relation between the point group lattice symmetry 
and the valley phase. As the parameter r is increased, another possibility is a topological phase 
transition to the X phase which does not respect the fourfold symmetry and thus distinguishes 
the square from the rectangular lattice. This is reflected in the response of the system upon the 
insertion of dislocations. In addition to the usual zero-energy modes resulting from a 7r-flux 
vortex, the spectrum shows a pair of zero-modes resulting from dislocations only if the Burgers 
vector b has a component in the x direction (Fig. S8). 
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Figure S7: Probing the X — Y valley phase, by considering one spin component of each 
Kramers pair, (a) Spectrum in presence of a 7r-flux and anti-7r-flux in a periodic system, (b) 
Spectrum in case of a dislocation and anti-dislocation in a periodic system, (c) Real space 
localization of zero-mode resulting from a single dislocation with Burgers vector b = e y . (d) 
The phase of the edge state hybridized with the dislocation zero mode displayed in (c), shows 
that it is a mode from the Y point, k = (0, re). Similarly, the dislocation zero-mode displayed 
in (f) hybridizes with edge state shown in (e), the phases of which indicate that this zero-energy 
state corresponds with the X point, k = (n, 0). Color-code is as in Fig. [T] The lattices contain 
33 x 33 unit cells. 

Let us conclude by addressing the robustness of the zero-modes in the topological non-trivial 
regimes, which is related to the stability of the corresponding phase. To this end, we multiply 
all the model parameters by Gaussian variables of width w = 10% and also introduce Gaussian 
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Figure S8: Illustration of the zero-mode bound to a dislocation appearing in the X phase 
decoupled from the edge states using linear combinations (left panel). The right panel shows 
the edge state, the phases of which indicate that these edge states indeed correspond with the X 
point (color-code as in Fig. [TJ). The lattice contains 30 x 33 unit cells. 

random chemical disorder, while preserving TRS. Moreover, we couple the spin-reversed sub- 
blocks by adding a nearest neighbor Rashba spin-orbit coupling 



to the Hamiltonian ( |ST| ), which breaks the reflection symmetry z — > — z about the plane. The 
results are displayed in Fig. S9. We find that finite Rashba spin orbit coupling, but not large 
enough to close the topological gap, preserves the real space localization of the zero-modes. 
More importantly, we observe that the stability of the modes in the valley phase is significantly 
smaller than in the Z 2 non-trivial phases. We note that the latter phase has a smaller gap and 
hence would already require a weaker Rashba coupling to close the gap. However, one would 
expect the valley phase to be less stable independently of the gap size, as only the lattice sym- 
metry accounts for its protection. To check this assertion, we also considered the valley phase 
in presence of a stronger next nearest neighbor hopping, which leads to an increase of the gap 
size. This analysis then confirms that the zero-modes in the valley phase are substantially less 
robust than those in the TRS protected phases. 




(S9) 



28 



M/B 



0, R=1.5 dis 



0.05 

qO.03 
Q0.02 
0.01 




O 
Q 



-0.5 0.0 
M/B=6, z=0, R=2 dis 




0.07 
0.06 
V)0.05 
0.04 
Q0.03 
0.02 
0.01 




-0.5 0.0 0.5 1.0 



M/B=8.81, z=0.7, R=0.5 dis 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 

0.25 

e M/B=10,z=0.9, R=l.dis 
0.030 
0.025 
0.020 
g 0.015 
D 0.010 
0.005 
0.000 

-1.0-0.5 0.0 



1 M/B: 

0.04 
t£>0.03 
§0.02 

0.01 



= 10.4, z=0.8, R=1.5dis 




-0.5 0.0 0.5 




Figure S9: The averaged (over 130 disorder realizations) density of states (DOS) for various 
values of M/B, z = B/B and R = R / B. Panels (a)-(b) and (d)-(e) correspond to the square 
lattice model, whereas panel (c) corresponds to the rectangular case. The midgap modes result- 
ing from dislocations are present even for strong Rashba coupling and disorder, as demonstrated 
by (a) R = 1.5 in the M phase, (b) R = 0.5 in the X - Y valley phase and (c) R = 1.5 in 
the X-phase. Sufficiently strong Rashba coupling closes the topological gap, as shown for (d) 
R = 2 in the M phase (a similar result is obtained for the X -phase) and (e) R = 1 in the X — Y 
valley phase. Even for comparable values of the bulk gap, the valley phase is substantially less 
robust to the Rashba coupling than the TRS protected M and X phases. 
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